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1 Introduction: Diquarks 

This is a status report of our work on quark/diquark effects inside compact astro- 
physical objects. It goes somewhat in excess of the results shown at the workshop. 

The word diquark is due to Gell-Mann in 1964 The idea of a two-quark corre- 
lation has now spread to many areas of particle physics, motivated by phenomenology, 
lattice calculations, QCD or instanton theory. Now there are some 1200 papers on 
diquarks, some of which were covered by a review in 1993 0. The speaker (S.F.) has 
an updated database with diquark papers. 

Although there is no consensus about diquark details, it seems certain that a (ud) 
quark pair experiences some attraction when in total spin-0 and colour-3*. Nuclear 
matter should also be subject to such pairing, maybe in some new way, e.g., as in 
the nuclear EMC effect |^. Quark pairing should also affect a (hypothetical) quark- 
gluon plasma (QGP), and by now there are more than 300 papers built on this idea. 
Many of these assume that the correlation is like the one between electrons in a 
superconductor, with diquarks being the Cooper pairs of QCD [Q. The notion of 
"superconducting quark matter" is due to Barrois in 1977 p. Current efforts owe 
much to the review paper by Bailin and Love ||^, and to work by Shuryak, Wilczek 
and collaborators 0, P). Diquarks in a QGP have also been analysed "classically", 
with thermodynamics, or field theory |T0|. 

Astrophysical diquarks gained popularity about a decade ago, when they were 



suggested to influence the supernova collapse and "bounce-off" [|r^, |T2|, |T3|, |T^, and 
to enhance the neutrino cooling of quark-stars. The latter effect is now subject to 
much research [T^ . 

Here we study some features of compact objects that would be sensitive to diquark 
condensation in a QGP. The form factor of the diquark correlation and the quark 
isospin (a)symmetry due to presence of electrons will be given special attention. 
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There are many situations where a QGP with diquarks might play a role: (i) In 
a quark star, which might appear as dark matter |T^, [18|; (ii) In a hybrid/neutron 
star, surrounded by a hadronic crust; (iii) In a supernova, or a 'hypernova' gamma- 
ray burster, where diquarks might trigger neutrino bursts and the bounce-off; (iv) 
In the primordial plasma at the Big Bang, where diquarks might have delayed the 
hadronisation. 



2 Formalism and Results 



We use the BCS theory of colour superconductivity @, 0, H- The gap A can be seen 
as the gain in energy due to the diquark correlation. Another gap, 0, is related to 
the quark-antiquark condensate. The thermodynamical grand canonical potential, 
Q, is minimised in its variables, resulting in an equation of state (EOS) and other 



relations. We follow the approach of [|1^ for the f2 and generalise it for isospin 
asymmetry between u and d quarks pO| 
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where the subtraction C = — fi((0g'^'^)^, 0; 0, 0, 0; 0) has been introduced to make (|l|) 
finite and to assure that pressure and energy density of the vacuum at T = /i = 
vanish. Thus at the boundary of a compact quark matter object, where the quark- 
condensate 00 changes, a pressure difference arises, which is necessary for confining 
the system, at least for small masses. 

Instead of the chemical potentials fiu, of u and d quarks, one uses ||2^ those of 
baryon number, = (/i^ + /id)/2, and isospin asymmetry, /i/ = (yU^ — /id)/2. Then, 
if beta equilibrium with electrons holds, /Xe = — 2/i/. The particle densities are given 
by riB = TLu + rid = —dQ/dfiB, nj = — na = —dQ/dfij and rig = —dQ/dfie = 
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Fi gure 1: The gaps 0, A and the isospin asymmetry chemical potential /i/ as functions of fiB for 
(a) Gaussian, and (b) NJL form factors. Solid lines are with isospin symmetry (and no electrons), 
Dashed lines are with beta equilibrium. Dotted lines in (a) are with isospin symmetry for Aq = 0.8 
GeV and t^™'^ = 0.33 GeV. Dashed double-dotted lines in (a) are with isospin symmetry for Aq — 0.8 
GeV and (Z)™"^ = 0.45 GeV (giving the same vacuum pressure as with the solid lines). 



— /ig/(37r^). Charge neutrality means ub + Srij — Qn^ = 0, and isospin symmetry 
ni = 0. 

The dispersion relations are 



E^ = Jq^ + {m + F^{q)(l)y and E± = {E^±fiB). 



where F{q) is a form factor for two-quark correlations. 

We use T = and m = 0, and start by finding the gaps (po and Aq that minimise 
Q. These values give the EoS as 

n(0o, Ao; /ifi, /i/, iJe;T = 0) = e - j^snB - fJ^inj - /ieWe = -P, (3) 

where e is the energy density and P the pressure. We use the form factors suggested by 
Schmidt et al. who fitted those to masses and decay properties of light mesons. 
The following special cases are considered: (i) Gaussian: = exp(— g^/A^); (ii) 

Lorentzian: F\q) = 1/[1 + (^/Al)^]; (in) cutoff NJL: F^{q) = 6(1 - q/ANJh)- Here 



bo = 0.33 GeV, A^ = 1.025 GeV, A^ = 0.8937 GeV and Anjl = 0.9 GeV [21 



We contrast this to A^ = 0.8 GeV used in Fig. 1 shows our results for the 

Gaussian and NJL form factors, and for two extreme cases: (i) isospin symmetry 
(equal numbers of u and d quarks), and (ii) beta equilibrium and charge neutrality. 
In the former case charge is balanced outside the object. In the latter, the electron 
fraction is small, and Ud ~ 2?t,„. Here ni is given in the lower part of the figures. We 
also test the sensitivity on the parameter A^. There is now the choice to either keep 
fixed, or to change it so that the vacuum pressure (0o'^'^)^/(4Gi) is kept fixed. 
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Figure 2: a) The mass ratio M/Msun vs. the radius i? of a compact object according the the TOV 
equation. The curves correspond to the same situations as in Fig. 1. 



The EoS is used as input to the standard Tolman-Oppenheimer-Volkoff equations 
22| for an equihbrium (spherical) fluid (see also p3|): 



dP{r) _ [e(r) + P(r)][m(r) + 47rr^P(r)] 
dr r[r — 2m(r)] ' 

m(r) = Att e{r')r'^dr' . (5) 







The equations are iterated in r, starting with some value eo at r = 0. The 
procedure stops when P = 0, which defines the radius R of the object, and M = m{R) 
is plotted vs. R, see Fig. 2. Each eo value results in one point in the plot, with unique 
values also of //^ and /i/. The graph becomes a backbending spiral, with eo and fis 
increasing along the curve. Solutions on left-going parts of the curve are unstable. 



3 Conclusions and Outlook 

The diquark form factor and the isospin (a) symmetry are important inside a quark/ 
diquark star. They can make a difference of up to two solar masses, as for its maximal 
mass. 

We now have results also for < T < A, where condensate-free states can occur 



in intervals of /i^ |]2J]- These correspond to an interaction- free plasma, and appear 
as "bumps" in the TOV curves. They might be relevant as phase transitions inside 
a collapsing system (supernova). The QGP might suddenly enter a diquark phase. 



whereafter it again becomes a free-quark state just before the bouncing We 
also intend to study such events for the Big Bang plasma. Maybe a transition into 
a diquark phase resulted in a release of neutrinos, just as with the cosmic microwave 
background. The neutrino energy would have been of the order of A, but now cooled 
to keV energies. 
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Discussion 



A. Thampan: Should one not see a signature (say in terms of phase transitions or 
so) in the P vs. p (EoS) relationships? 

Fredriksson: For some extreme parameter values bumps occur in, e.g., M vs. R 
and P vs. e, and certainly also in a P vs. p, although we did plot the latter. For 
the case of beta equilibrium, these transitions are smoother in the TOV graphs. The 
bumps are clearer for T > 0, as mentioned above. 

J.E. Horvath: Do diquarks disappear suddenly in this model? Could you identify 
what makes the matter self-bound? 

Fredriksson: Not suddenly in time, because these are equilibrium solutions. But 
diquarks disappear "suddenly" when other parameters change, e.g., (see "NJL"). 
A strange star is self-bound by gravity at high masses, but also due to our "vacuum 
pressure" (C). 

Unknown: Are the stars beyond the mass peak stable? 

Fredriksson: Our approach does not tell, but states with ^ < are usually 
prescribed to be unstable [|^ . 
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